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Abstract 
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O 

O 

t^^ ' We calculate the low-momentum N-N effective potential obtained in the 

O' 

^r I OBE approximation, inside a nuclear plasma at finite temperature, as de- 

scribed by the relativistic a-uj model. We analyze the screening effects on 
O I the attractive part of the potential in the intermediate range as density or 



c 



temperature increase. In the long range the potential shows Friedel-like os- 
cillations instead of the usual exponential damping. These oscillations arise 



X 

C^ I from the sharp edge of the Fermi surface and should be encountered in any 

realistic model of nuclear matter. 
21.65. +f,21.30.+y,21.60.Jz 



Typeset using REVTgX 



In the usual approaches to the study of nuclear matter near saturation density, the phe- 
nomenological potentials are basic ingredients for the description of the nuclear interaction 
in a Schrodinger-dynamical framework ||T| . Potentials obtained from the exchange of single 
pseudoscalar, scalar, pseudovector and vector mesons with different theoretical approaches 
lead to satisfactory quantitative predictions of the observed properties of two-nucleon sys- 
tems. A simplified description of the main features of the N-N interaction can be obtained 
from the exchange of scalar and vector mesons only [^ . 

For higher densities (owing to the Pauli principle) relativistic effects become essential, 
not only in the description of the interaction itself, but also for the analysis of the particle dy- 
namics. In this case, the Lagrangian approach becomes the natural framework for the study 
of the nuclear plasma |^. The solution of appropriate Lagrangian models in the relativis- 
tic Hartree approximation (RHA) provides a satisfactory picture of the thermodynamical 
behaviour of relativistic nuclear matter at finite temperature |^ . Moreover, the relativistic 
meson propagators in vacuum obtained from this approach allow for the calculation of the 
N-N interaction potentials |^. 

Inside nuclear matter, the polarization effects introduce important changes in the form 
of the relativistic meson propagators. Consequently, the N-N interaction potentials inside 
the plasma are also strongly modified by the screening. Their behaviour as a function of the 
thermodynamical state provides a very explanatory visualization of these effects and suggest 
the existence of new collective phenomena. 

In this letter we report some results of a study of the N-N interaction potentials ob- 
tained in the one-boson exchange model (OBE) inside symmetric nuclear matter at finite 
temperature, and analyze some phenomena related to the screening. 

The Lagrangian model used for this calculation describes the nuclear interaction in terms 
of scalar-o" and vector-u; mesons exchanges . Although such a simplified model is not able 
to account for the whole richness of actual nuclear matter, it gives (when solved in RHA 
0) an acceptable description for its thermodynamical behaviour. Moreover, as mentioned 
above, the involved mesons reproduce the main qualitative features of the nuclear interac- 



tion in vacuum. In fact the fictitious cr-meson, which provides the attractive part of the 
potential at intermediate range, is introduced as a simple parameterization of the correlated 
2 TT-exchange contribution to the N-N interaction. It has not been proven that such a pa- 
rameterization is also possible inside the nuclear medium. Nevertheless, the behaviour of the 
screened potential with the thermodynamical state found here is qualitatively similar, in the 
intermediate range, to the one found in a 2 7r-exchange calculation P]. Moreover, as we shall 
see, the medium effects on the long-range behaviour of this potential are mainly dominated 
by the singularities of the matter polarization contributions, and are rather independent of 
the details of the basic interaction. They should be present in more realistic analysis which 
take into account all the relevant meson exchanges. 

When solving the model in the RHA, the values of the constants in the Lagrangian 
are fixed as follows: The meson and fermion masses are fixed to their "physical" values 
fia = 550 MeV, /io; = 783 MeV and m = 939 MeV. For the coupling constants, we 
choose the values which lead to a satisfactory fit of the saturation properties in the RHA: 
gl = 183.3(/io-/?7i)^ and (^5 = llA.7{^t^/my. From these values, saturation is attained at a 
Fermi momentum of the nucleon PfQ = lA2fm~^, with a binding energy Ef, = —15A6MeV. 

In going beyond RHA, the analysis of the small perturbations of the fermion distribu- 
tion and meson fields around the Hartree equilibrium gives the expressions for the meson 
propagation equations inside the plasma, which can can be written in a matrix form as 
D{k).J]i{k) = 0, where Si is the column matrix of the components of the scalar and vector 
perturbing fields, and D{k) is a 5x5-matrix containing the scalar (Hq.), mixing (H^^^) and 
vector (H^'^) polarization tensors, which are functions of the thermodynamical state and 
include the renormalized vacuum contributions 0. At T=0 they reduce to the one- loop 
meson polarizations ||10[. The propagator matrix for the mixed scalar- vector field inside 
matter is given by P,pJ]] G{k) = — D{k)^^. At zero density and temperature, the meson 
fields decouple from each other, and the components of this matrix reduce to the one-loop 
scalar and vector propagators in vacuum. 

The first step in the calculation of the screened two-nucleon potential is the derivation 



of the relativistic one-boson {a + u) exchange amphtude diagrams from the propagator. 
Now, the difference with the calculation in vacuum is that the meson propagators are 
dressed by the medium, and mixing between both meson fields appears |p. Therefore, the 
usual Feynman rules are slightly modified in the present case: We must introduce a iTa 
factor at each vertex (Fq is the coupling matrix defined as F^ = 'y^g^^ for a = /i from 
to 3, and F4 = g^ ), and a dressed boson propagator matrix iG{k) for each internal 
boson line. Moreover, in the present calculation we are interested in the structure of the 
propagator matrix on the k^ = axis only, where there are two poles associated to 
the "tachyonic" branches, coming from the vacuum polarization terms |^. Such poles are 
spurious because they arise at large values of q, where the point-particle approach fails, and 
the nucleon structure should be taken into account. This is done through the introduction 
of phenomenological monopolar form factors: 

fa{k) = {Al - ^l)l{Al - e) (1) 

(a = cr,u;); 

at each vertex of the loop and boson-exchange diagrams. This amounts to multiplying 
each squared coupling constant by the corresponding form factor. With this prescription the 
spurious "tachyonic" branches (and the associated poles in the propagator matrix) disappear 

i- 

Under these conditions, the relativistic OBE amplitude takes the form 

■ {zG{k)r^ { [ xl' u{p'2, 4) ] (*r)„ [ u{p2, S2) X2 ] } (2) 

where pi and p2 are the four-momenta for the incoming quasi-nucleon states and p'l 
and P2 correspond to the outgoing quasi-nucleon states, whereas k = pi — p[ = p'2 — P2 is 
the transferred four-momentum. The interacting quasi- nucleons have effective mass M (as 
given by the RHA approximation) and spins 5*1 = |(Ti and S2 = \^2 ■ The indices m and n 
run from to 4 . Finally u{p, s) and x ^ire the Dirac cuadri-spinor and the isospin wave 
function, respectively. 



The potential (OBEP) is obtained by eliminating the wave functions of the initial and 
final states in the amplitude Eq. (|2|) taken in the center-of-mass system. From this ex- 
pression we obtain the non-relativistic potential by performing an expansion in the nucleon 
momenta and keeping only the second-order terms. We have corrected the OBEP according 
to the Blankenblecker-Sugar prescription |T^, which includes the requirement of "minimal 



relativity" g as defined by V{q',q) = {M/Eg>y/^VoBEp{q',q) {M/Egf'^ where q and q' 
are the CMS initial and final momenta of the nucleons, respectively, Eq = \JM'^ + q^^ and 



Eqi = Y M^ -|- q'^. Moreover, the above prescription implies taking the static (A;° = 0) 
limit. After Fourier transformation one obtains the potential in coordinate space: 

V(f) = K(r) - ^( VV2(f) + ^(r) V^) + VLsir) L-S + Vss{r) a, ■ a^ + Vrir) Su (3) 

where L = f Ap is the orbital kinetic momentum, S = (1/2)((Ti + 0^2) is the total spin 
operator and 5*12 = ^(f?i • x) ((T2 ■ x) — ( ai • ^2 ) is the tensor operator. 

The second term in Eq.(0) is a non-local component which gives a small contribution 
to the potential. We shall omit here the study of this component. The other terms (cen- 
tral, spin-spin, spin-orbit and tensor components) are now functions of the interparticle 
distance, as well as of the thermodynamical state of the plasma (density and temperature). 
We must emphasize that the non-relativistic limit concerns the dynamics of the two inter- 
acting nucleons and the neglecting of retardation effects. However, the description of the 
thermodynamical state of the plasma in the RHA remains fully relativistic. Also, no low-q 
approximation has been done for the meson propagators. 

In the analysis of the potential, the values of the coupling constants and cut-off pa- 
rameters in the form factors have been fixed in order to fit the deuteron and low-energy 
phenomenology data 0: 

gl = 8.7171 9^ = 25 A^ = 2.0GeV A^ = lAGeV (4) 

However, in calculating the underlying thermodynamical state, the coupling constants 
have been fixed to the above mentioned values which fit saturation in RHA. 



Figure (1) is a picture of the central component of the potential at T = 0, at saturation 
density (solid line) and 2.4 times these density (dashed-dotted line), for symmetric nuclear 
matter [Pj is the Fermi momentum of the nucleons). The same component, calculated at 
zero density, with (Vpoi) and without [Vyac) vacuum polarization contributions has also been 
plotted. The very short-range region (r < 0.6/m) must be discarded because the potential 
behaviour is dominated there by large-q values, and the non-relativistic (low-q) limit breaks 
down ]ll|. Beyond this distance, a hard-core and a potential well appear both in vacuum and 
at finite density. At zero density, the vacuum polarization effects enhance the slope of the 
hard-core and increase the depth of the potential well by an important amount. The slope 
of the repulsive-vector Yukawa component is strongly raised by the vacuum polarization, 
whereas the attractive-scalar component remains nearly unaffected by these effects. This 
explains the observed behaviour in vacuum. The matter polarization effects increase the 
range of the vector component, and therefore the depth of the well is reduced as density 
grows. It disappears shortly above the saturation density |11|. Similar results are obtained 



if the attractive a contribution to the N-N potential is replaced by the 2 vr-exchange [^]. 

At larger distances, the screening effects of the medium introduce important qualitative 
new features. Whereas in vacuum the potential shows an exponential damping with the 
distance, an oscillatory behaviour appears at finite density, whose amplitude is damped as 
an integer power of the distance. These are Friedel oscillations, similar to those encountered 
in many low-temperature Fermi systems [|13[. Mathematically, these oscillations arise from 



the fact that the matter polarization contributions to the screened meson propagators show 
singularities in their derivatives at g = 2Pf (Kohn singularities [0). After Fourier trans- 
formation, such singularities introduce oscillations in the r-space potential. More detailed 



analytical calculations |T^, [|1^] show that, in this region, the potential can be decom- 
posed into a and u Yukawa-like components (which dominate at shorter distances), and 
this Friedel-like component, which determines the large-r features and is long-ranged and 
oscillatory. From a physical point of view, the Kohn singularity is associated to the sharp 
character of the Fermi surface at T = and is rather independent on the details of the 



interaction. Consequently, the long-range oscillations of the screened interparticle potential 
must be an universal feature of interacting degenerate Fermi systems. Indeed, this has been 
found for degenerate non-relativistic [0 and relativistic |T^ electromagnetic plasmas, for a 



QCD plasma |T8[, and for the screened one-pion exchange potential in a relativistic nuclear 
plasma []T3 . 



In the calculation of Ref. f^, where the 2 vr-exchange is explicitly considered, no Friedel 
oscillations have been found. In fact, this paper was concerned with the intermediate range 
of the screened N-N potential, and the complete effects of the matter polarizations were 
not included in the meson propagators. Such effects go beyond the effective meson mass 
variation with density, considered by the authors. Therefore, it is not surprising that they 
do not obtain oscillations in the long range. 

As a test of the small sensitivity of the Friedel behaviour to the details of the interaction, 
we have plotted in Fig.(l) the central component of the potential at saturation density and 
T = (dotted line), now using for the potential the values of the coupling constants which 
fit saturation in the RHA, which are very different from the ones used in the previous 
calculation for the same thermodynamical state. We observe important differences in the 
intermediate region in both cases but, in the long range, the oscillations of the potential are 
only slightly affected by the changes in the intensity of the couplings. 

In Fig. (2) the spin-spin component of the potential is plotted at T = 0, in vacuum 
{Pf = 0), at saturation density {Pf/m = 0.3), and 2.4 times saturation density {Pf/m = 
0.4). A glance to this drawing shows also a Friedel-like oscillatory behaviour whose amplitude 
increases with density. (The spin-orbit and tensor components of the potential show also 
a similar behaviour: in all cases the amplitude of the oscillations increase as density grows 
pl|). At saturation density, the most important effect in the long-range appears on the 
central potential component, where the amplitude of the first oscillation reaches 2MeV 
for a distance of around 2fm. Nevertheless, for higher densities the oscillations of the 
spin-spin component reach comparable maxima. Indeed, for Pf/m = 0.4 the amplitude of 
the first oscillation in the central component is ~ 8MeV at Vmax ~ 1.5/m (see Fig.(l)), 



whereas the spin-spin component reaches lOMeV at the same distance (see Fig. (2)). The 
position of the first maximum of the oscillations, r^ax , is related to the Fermi momentum 
through |T^ Ppf-Tmax ~ TT. The ratio between the position of the first maximum of the 
oscillation and the mean interparticle distance d = {37r'^/2P?Y^^ is nearly constant with 
density {rmax/d ^ 1.28). 

At finite temperature, the slope of the hard core is slightly reduced, as showed in Fig. (3), 
where the central component has been plotted at saturation density and various tempera- 
tures. In the intermediate range, the effects of temperature on the central potential reduce 
the depth of the well. The well disappears for temperatures beyond 40 MeV. This can be 
interpreted in terms of the modifications introduced on the Yukawa-like components by the 
temperature in this region. As temperature increases, the ranges and intensities of these 
Yukawa components are modified in such a way that in the balance, the vector-repulsive 
part becomes dominant and the well disappears, (see Refs. [|15| and |1^ for a more detailed 
analytical study). 

Concerning the long-range behaviour of the potential, as temperature increases the Fermi 
distribution function becomes smooth and the Kohn singularity (and the associated oscil- 
lations) disappears in all the components of the potential. This is indeed observed in Fig. 
(3). The amplitudes of the oscillations decrease with temperature. As can be checked an- 
aljd;ically [1I5| , |T9| this amplitude at a fixed distance and constant density is exponentially 
damped with temperature. The oscillations in all the components disappear for tempera- 
tures between 40 — 80MeV. Beyond this, the potential becomes exponentially damped with 
the distance. 
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FIGURES 
FIG. 1. Central component of the potential at T = at zero density with (Vpoi) and without 

(Vvac) vacuum polarization contributions at saturation density (Pf/m = 0.3) and 2.4 times 

saturation density (Pf/m = 0.4). The dotted line (with Pf/m = 0.3) corresponds to the 

central component of the potential at saturation density obtained from different values of the 

model parameters. 

FIG. 2. Spin-spin component of the potential at T = for the same values of the Fermi 
momentum as in Fig.l. 

FIG. 3. Central component of the potential at saturation density for different temperatures. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9407002vl 



This figure "fig2-l.png" is available in "png" format from: 
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